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Abstract 

We present a family of time-dependent solutions to 2+1 gravity 
with negative cosmological constant and a massless scalar field as 
source. These solutions are continuously self-similar near the central 
singularity. We analyze linear perturbations of these solutions, and dis¬ 
cuss the subtle question of boundary conditions. We find two growing 
modes, one of which corresponds to the linearization of static singular 
solutions, while the other describes black hole formation. 


1 Introduction 

In recent numerical simulations jlj of gravitational collapse of a massless 
scalar field in 2+1 dimensions with negative cosmological constant A, thresh¬ 
old solutions were observed, exhibiting features characteristic of critical col¬ 
lapse m, namely power-law scaling and continuous self-similarity (CSS). It 
has been argued by Garfinkle p[] that this behavior is governed by an exact 
CSS solution of the A = 0 theory. To show that such a solution is indeed 
a critical solution, one should investigate near-critical configurations to see 
whether they describe black hole formation. However, as a strictly negative 
cosmological constant is necessary for the existence of black holes in 2+1 
dimensions Q, one should first extend the Garfinkle solutions to A < 0, a 
hard problem which has not yet been solved in generality. 

In this talk, we shall present a new class of A = 0 CSS solutions which 
can be extended to A < 0, and show that a class of nearby solutions do 
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indeed describe black hole formation. After reviewing the Garfinkle solu¬ 
tions, we describe the new CSS solutions and their extension to A / 0 
quasi-CSS solutions. We then perform the linear perturbation analysis in 
the background of these threshold solutions. Imposing appropriate bound¬ 
ary conditions, we obtain discrete modes describing black hole formation, 
and discuss the determination of the critical exponent. 


2 CSS and quasi-CSS solutions 


The Einstein-massless scalar held equations with cosmological constant are 


G, w - A = d^(f)d v (p - (1/2)^0 a </>5a</), V 2 0 = 0. (1) 


Assuming rotational symmetry, we make the double null ansatz: 

ds 2 = e 2a{u ’ v) dudv - r 2 {u, v)dd 2 , </> = v). (2) 


The held equations then take the form 

r,uv = (A/2 )re 2a , (3) 

2a, uv = (A/2)e 2fr - <t>, u (j>,v, (4) 

2 v, u r,u ~ r )U u = r4 > 2 u , (5) 

2o>J> - r >vv = rcj) 2 v , (6) 

2 r<j) >uv + r >u (f> )V + r >v <j> >u = 0 . (7) 


Assuming A = 0, Garhnkle found Q the following exact CSS solutions 


ds 2 

(t> 


-A 



dudv — (v + u) 2 d9 2 , 


— 2cln(Vi + \[— 


( 8 ) 

(9) 


depending on an arbitrary constant c and a scale A > 0. These solutions are 
continuously self-similar with homothetic vector (uda + vdy) and, for c ~ 1, 
agree well near the singularity with numerical results of Q. However, it has 
not been possible up to now to extend them to solutions of the full A < 0 
equations with both CSS behaviour near the singularity, and the correct 
AdS behaviour at spatial infinity. 

From the Garhnkle class of CSS solutions (D-® , we derive a new class 
of CSS solutions by an infinite boost it —> D —> Ah (A —> 0), arriving 

at 

ds 2 = — A(—u / v) c2 dudv — u 2 d0 2 . (10) 
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These new solutions obviously inherit the CSS behaviour of the original 
Garfinkle solutions. The transformation to new null coordinates u = — (— u) l+c , 
v = fr c leads to the simpler form of these solutions 

ds 2 = dudv — (—u) 2a dd 2 , ^ = — caln(— u) ( 11 ) 

(a = 1/(1 + c 2 )), which may be checked to solve exactly Eqs. (!)-©• 

The generic solution ([TT| ) has four Killing vectors. While this spacetime 
is devoid of scalar curvature singularities, the study of geodesic motion shows 
that for c 2 < 1 only radial geodesics can be continued through the null line 
u = 0 , while nonradial geodesics terminate at the singular point u = 0 , v —> 
Too. On the other hand, for c 2 > 1 nonradial geodesics terminate on the 
null line u = 0. So the c 2 = 1 solution is an extreme solution separating the 
c 2 < 1 solutions with a point singularity from the c 2 > 1 solutions with a 
null line singularity. 

Now we proceed to extend the second class of A = 0 CSS solutions (ill:) 
to exact quasi-CSS solutions of the full A < 0 equations (the self-similarity 
being then broken by the cosmological constant). We make the ansatz 


ds 2 = e 2v ^dudv — (— u) 2a p 2 {x)d6 2 , cj) = — ca ln(— u) + i p(x) (x = uv). 

( 12 ) 

Inserting this ansatz into the field equations (!|)-(@) leads to the system 


xp" T (1 T a)p' 

= (mpe 2 \ 

(13) 

2(xu' / T v') T i/j'ixip' — ca) 

= (A/ 2 )e 2 -, 

(14) 

x 2 [—p' T 2pV — pip' 2 ) T 2 ax(—p' + p{v' T cip')) 

= 0 

(15) 

-p" + 2pv - p r df 2 

= 0 

(16) 

2 x(pi/;'y T (2 T OL)pi\)' 

= cap'. 

(17) 

A simple first integral is 

p = e u+c ^. 


(18) 


The numerical solution of this system with the boundary conditions 


p(0) = l, v(0) = 0 , 4>{0) = 0 (19) 

leads [||] to a unique extension which reduces to (|ll]) near the singularity 
u = 0 and which, for x —> x\ < 0, is asymptotic to the AdS 3 (anti-de-Sitter) 
spacetime 

ds 2 = (X 2 T 1 )dT 2 - i - X 2 d0 2 . (20) 
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3 Perturbations 


To ascertain whether the quasi-CSS solutions of the preceding section are 
threshold solutions for black hole formation, we study their linear pertur¬ 
bations, along the lines of the analysis of [| 6 |. The relevant time parameter 
in critical collapse being the retarded time |] T = — ln(— U) = — aln(— u), 
we expand these perturbations in modes proportional to e kl = (—u)~ ka . 
Keeping only one mode, we decompose the perturbed fields as 

r = (— u) a [p{x) + A(— u)~ ka r(x)\, ( 21 ) 

a = u(x) +X(-u)~ ka a(x), ( 22 ) 

4 = —caln(—u) + ip(x) + \(—u)~ ka ^>(x). (23) 

The linearization of the Einstein equations (|j)-(|7j) in the small parameter 

A leads to the fourth order system for the perturbations f, a, 4>: 

xr" + (1 + a — ka)f' = (A/2) e 2iy (f + 2pa), (24) 

2 xa" + 2(1 — ka)a = Ae 2u a — ( 2x^ ' — ca)4> + kai/j'cf), (25) 

— (1 — k)xr' + ((1 — k)xv 1 + (k/2){2a — 1 — ka))r + pxa' — k(xp' + ap)a = 
—p{cx(p — k(ca — xip')(f)) — cxi/j'r, (26) 

2 (p’a' + v'r') - f" = il}'(2p& + ip'r), (27) 

2 xp4>" + {2xp + (2 + a — 2kap)(\) — kap'cj) + ( 2x — ca)r' 

+{2xi\}" + (2 + a - ka)^')f = 0 . (28) 

This system admits the exact, spurious solutions (gauge modes) 

= (—x) p +p', = (— x) v+ v' — ^K(— x) p+ ~ l , 4^ = (— x) p+ i/j', (29) 

= xp’(x) + ap , a^ ] = xu 1 + (ffi = xij/ - ca , (30) 

(p± = 1 ± ka) generated respectively by the small gauge transformations 
v —> v — Xv p+ and u —> u — X(—u) p ~. So, up to gauge transformations, the 
general solution of this system depends only on two integration constants, 
which should be determined by enforcing appropriate boundary conditions. 

The three boundaries of our quasi-CSS solutions are u = 0, v = 0 (or 
x = 0), and x = x\ (the AdS boundary). On the singular boundary u = 0, 
it seems natural to impose that the perturbed metric component r does 
not diverge too quickly, which would conflict with the linear approximation. 
However we shall argue that it may be preferable to replace this condition 
(which is subject to some ambiguity) with a boundary condition on the 
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apparent horizon. On the second boundary v = 0, we shall require that the 
perturbed solution matches smoothly the original quasi-CSS solution 

r(0) = 0, <x(0) = 0, 0(0) = 0, (31) 

and moreover that the perturbations be (as the unperturbed quasi-CSS so¬ 
lution) analytic in v. This last condition ensures regularity on the line v = 0. 
Finally, the perturbations should grow slowly (i.e. not faster than the orig¬ 
inal fields) as the AdS boundary x = x\ is approached. It turns out that 
this last condition is always satisfied up to gauge transformations. This is 
due to the fact that the perturbed scalar field becomes negligible at spa¬ 
tial infinity, so that the asymptotic behavior is that of AdS 3 (the effect of 
gauge transformations is to shift the value x± for which the AdS boundary 
is attained). 

To enforce our boundary conditions on the line v = 0, it is enough to 
know the behavior of the solutions of the linearized system near x = 0. Up 
to gauge transformations, this is given in terms of two integration constants 
B,C by 

f{x) ~ B(-x) + AC(—x) 1+(fc ~ 1/2) ", (32) 

a(x) ~ —4>(x) ~ —(1 + a — ka)A~ 1 B + 0(x) 

-(1 + a/2)(l + (k — 1/2 )a)C(-x) (k ~ 1/2)a . (33) 

for k 7 ^ 1/2 (in the degenerate case k = 1 / 2 , the power law (—aj)( fc-1 / 2 ) a is 
replaced by a logarithm). From these behaviors, we find that the smooth 
matching conditions at v = 0 are satisfied if either 

(7 = 0, k = c 2 + 2, (34) 

(in this case the analyticity condition is also satisfied), or B = 0 which, 
taking into account the analyticity condition, implies 

5 = 0, k = n(c 2 + 1) + 1 / 2 , (35) 

where n is a positive integer. In this last case the boundary condition that 
the perturbed r does not diverge too fast (that is, no faster than in the BTZ 
case, see Sect. 4) when u —>■ 0 with x fixed selects the eigenvalue n = 1, i.e. 
k = c 2 + 3/2. 

4 Black hole formation 

Black hole formation is characterized by the appearance of a central singu¬ 
larity hidden by an apparent horizon. The definition of these notions can be 
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Figure 1: Perturbative singularity (thick curve) and apparent horizon for 
solutions (|34|) (left) and (|3|) (right) 


quite ambiguous in linearized perturbation theory. For our present purpose 
we shall identify the central singularity with the coordinate singularity of 
the linearized perturbed metric 

y/g = e 2a r = (—u) a e 2u ^ x ' > [p(x) + (—u)~ ka (r(x) + 2 p(x)d(x))\ = 0, (36) 


and define the apparent horizon by 

r> = (-u) 1+ V + {-u)- ka r') = 0. (37) 

In the case of the solution (|34j), for A B < 0, a spacelike coordinate 
singularity appears at v = ( 1 . hidden behind a preexisting spacelike apparent 
horizon, which seems to be eternal (Fig. 1). Such a situation does not 
seem to correspond to collapse, but rather to a static configuration. Indeed, 
consider the exact static solution of the Einstein equations (]3])-([7j) with A = 0 

B 

ds 2 = Ar c2 (dt 2 — dr 2 ) — r 2 d6 2 , 0 = clnr. (38) 

Putting r = (—u) a + Xv, t = — (— u) a + Xv, and expanding in powers of the 
small parameter A, we get to first order 

ds 2 0 = A(l — Xc 2 {—u)~ a ~ l )dudv — (—u) 2q (1 — Xc 2 {—u)~ a ~ 1 ) 2 d9 2 
4> = — cln(— u) + Ac(— u)~ a ~ 1 x, (40) 


which agrees with the first order expansions (21)-(p3|) for A = 0 and ka = 
a + 1, i.e. k = c 2 + 2 . We conjecture that the full A < 0 perturbed solution 
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with k = c 2 + 2 results from a similar expansion of the “non-topological 
soliton” singular static solutions (generalizing (^|) to A < 0) || 


ds 2 


A\p- p + \ l / 2+a \p- p_\ l / 2 ~ a dt 2 
■ dp 2 

+ 4A(p-p + )(p-p-)’ 9 


A 


Ip-p+l^lp-p-l^de 2 


1-4 a 2 


In 


/ |p- p+h 

vip-p-|/ 


(41) 


(with 2a = (c 2 — 2)/(c 2 + 2)). This conjecture is certainly true in the 
sourceless case c = 0, for which these solutions reduce to the BTZ black 
holes 

dr 2 

ds 2 = ( r 2 /l 2 - M)dt 2 - (r , 2//2 _ M) - r 2 dd 2 , $ = 0 (42) 

(. 1 2 = —A -1 ). The BTZ metric can be rewritten in the double-null form j5] 

ds 2 = (x + uv/Al 2 ) 2 ^ iudv ~ U + ( M / 4 ) u ) 2d6)2 ]> ( 43 ) 


which corresponds to the first order expansion (21)-(23) for p = e v = (1 + 
uv/41 2 )- 1 , (/) = 0, A = M/4, k = 2 (the linearized solution is exact in this 


case). 

Consider now the second solution (fill). For C > 0, a coordinate sin¬ 
gularity and an apparent horizon, both spacelike, appear simultaneously at 
v = 0 , and converge towards the AdS boundary x = x\ for a finite value tah 
of the horizon radius (Fig. 1). The fact that this apparent horizon is born 
together with the singularity suggests that it is not automatically regular. 
Indeed, the curvature invariant 


i? + 6A = 4e 2,7 


( 44 ) 


evaluated on the apparent horizon, generically diverges on the original null 
singularity u = 0 because ^,v\ah x^ 1 (x —> 0). So the apparent horizon 
can be regular at birth only if 

= 0) = 0, (45) 

which leads to an eigenvalue equation relating n and c 2 . While there is some 
ambiguity in the linearization of (|4^) , it seems to lead jB| to a single solution 
n = 1 , c 2 ~ 1 . 

Near-critical collapse is characterized by a critical exponent 7 , defined 
by the scaling relation 

Q oc |p — p* | S7 (46) 
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for a quantity Q with dimension s depending on a parameter p (with p = p* 
for the critical solution). In the perturbative approach, the choice Q = rAH 
(the apparent horizon radius, s = 1) and p — p* = A (the perturbation 
amplitude) leads to 7 = 1/k ||. For the static BTZ black hole (solution 
(|34]) with c 2 = 0 ), we thus obtain 

1(BTZ ) = 1/2; (47) 

in agreement with previous analyses. In the case of genuine scalar field 
collapse (solution (|35|) with n = 1 ) we obtain for the preferred value c 2 = 1 
the critical exponent 

7 = 2/5. (48) 

This differs from the two conflicting values (7 ~ 1.2 and 7 ~ 0.8) found in 
the numerical analyses of [l|], but is of the same order of magnitude. 


5 Conclusion 

We have constructed a one-parameter (c 2 ) family of exact time-dependent 
solutions to the Einstein equations with a negative cosmological constant 
and a massless scalar field as source. These solutions are continuously self¬ 
similar near the singularity, and asymptotically AdS at spatial infinity. We 
have then analyzed the linear perturbations of these solutions, and discussed 
the subtle question of boundary conditions. These led to two possible grow¬ 
ing modes. The first mode, (j34|), is characterized by the presence of a 
presumably eternal apparent horizon. This is simply the linearization of the 
static BTZ black hole in the sourceless case, and very probably of static sin¬ 
gular solutions in the general coupled case. The second mode, (^), exhibits 
the characteristic features of black hole formation, i.e. the simultaneous ap¬ 
pearance of a spacelike singularity and its shielding apparent horizon, which 
then expands up to a finite size. The condition of regularity at birth of 
this apparent horizon selects a unique value of the parameter c 2 , leading to 
a value of the critical exponent which is of the order of the values derived 
from the numerical simulations. While these results are rather satisfactory, 
we think that further analytical work on this simple 2+1 dimensional model 
should further clarify the situation and help us to understand better the 
phenomenon of black hole formation. 

We would like to thank Marco Cavaglia for discussions. 
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